Decay properties of the one-particle Green function in real space and imaginary time 
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The decay properties of the one-particle Green function in real space and imaginary time are 
systematically studied for solids. I present an analytic solution for the homogeneous electron gas 
at finite and at zero temperature as well as asymptotic formulas for real metals and insulators that 
allow an analytic treatment in electronic-structure calculations based on a space-time representation. 
The generic dependence of the decay constants on known system parameters is used to compare 
the scaling of reciprocal-space algorithms for the GW approximation and the space-time method. 
Copyright 2000 by The American Physical Society. 
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I. INTRODUCTION 



The GW approximationEI for the electronic self-energy 
is known to produce quasiparticle band structures in very 
good agreement with experimental data for a wide range 
of crystalline Materials that include semiconductorsp'cl 
simple metalSjQ and even transition metals.cl It corre- 
sponds to a summation of ring diagrams to infinite or- 
der and may be thought of as an extended Hartree-Fock 
scheme with dynamically screened exchange. Despite 
the long-recognized success for bulk systems there are 
still comparably few applications to more complicated 
geometries, however, because of the rapid increase in 
computational cost. Ilaless additional simplifications like 
a plasmon-pole modcfl for the screened Coulomb inter- 
action are employed, conventional reciprocal-space algo- 
rithms scale like N'^/b'^ for a calculation of the complete 
spectral function, where iV is a measure of the system 
size, such as the number of atoms in a bulk material, 
and b is inversely aroportional to the number of fre- 
quency mesh points.tl The scaling problem is further com- 
pounded if diagrammatic vertex corrections beyond the 
GW approximation are included.Q The design of more 
efficient implementations is therefore a pressing task. 

In order to achieve a more favorable scaling, Rojas et 
at recently proposed an algorithm based on an alterna- 
tive representation flfnthe Green function in real space 
and imaginary timeE'EI This approach allows the calcula- 
tion of the irreducible polarizability and the self-energy as 
simple products rather than numerically expensive con- 
volutions. Furthermore, it directly exploits locality, i.e., 
the fact that the Green function decays to zero as its spa- 
tial or temporal arguments move apart. As a result, with 
a fixed cutoff energy this so-called space-time method 
only scales like N'^ /j^b, where 7 denotes the exponential 
decay constant in real space. Logarithmic terms are ne- 
glected. Under the assumption that the frequency mesh 
maps one-to-one onto the relevant region of the imagi- 
nary time axis, b can there be identified with the inverse 
period ^ksT in a finite-temperature formalism and with 
the exponential decay constant otherwise. The locality 
of the Green function is a maniijestation of the more gen- 
eral "nearsightedness" principleEj and due to destructive 



quantum interference in many-electron systems. 

While the behavior of the Green function in recipro- 
cal space and iiie frequency domain has been analyzed 
in great detail,tll there has so far been no comprehensive 
study of the decay properties in real space and, for zero 
temperature, on the imaginary time axis. In fact, even 
qualitative features are not always recognized correctly: 
Ref. 1^, for instance, describes the asymptotic tail of the 
Green function as being proportional to |r — r'|^^. While 
this assertion is correct for the homogeneous electron gas 
at zero temperature, I argue here that real materials in 
general exhibit exponential rather than algebraic decay. 
Attempts to model the self-energy rUf ips results obtained 
for the homogeneous electron gasJlJO which are moti- 
vated by seemingly uniiiersal features in the short-range 
part of the nonlocalitwij and the success of the locak 
density approximationliS in density-functional theory,EZI 
should appreciate such fundamental differences. 

The first objective of this paper is to systematically 
study the decay properties of the one-particle Green func- 
tion in solids, separately for metals and insulators, and to 
obtain analytic expressions for the asymptotic behavior 
in real space and imaginary time. Besides intrinsic inter- 
est, these can be directly exploited in electronic-structure 
calculations. In the space-time method, for instance, an 
analytic treatment of the asymptotic tail on the imagi- 
nary time axis, in combination with a customized Gauss- 
Legendre grid, reduces the number of mesh points where 
the Green function must be evaluated numerically by one 
order of magnitudeJl3 I show here that the exponential 
fit proposed in Ref. ^ is only valid for systems with a 
finite band gap, however, while metals at zero temper- 
ature instead require an algebraic fitting function. The 
second objective is to relate the decay constants 7 and 
b to known system parameters. This is then used for a 
more detailed scaling comparison between the space-time 
approach and conventional implementations. j— . 

Due to the central role of locality in 0{N) methodsE£l 
within density-functional theory, the spatial decay rate 
of the one-particle density matrix, which is a special el- 
ement of the-Green function, has recently received some 
attentionBj'EJ Physically, the nonlocality of the density 
matrix stems entirely from the delocalized quantum char- 
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acter of the wave functions, while the Green function also 
describes the actual transport of particles. The results 
obtained for the density matrix can therefore not be gen- 
eralized in a straightforward way, although they are re- 
covered here if the imaginary-time argument approaches 
zero from below. Furthermore, the one-particle density 
matrix contains no information about the behavior of the 
Green function on the imaginary time axis. 

This paper is organized as follows. Sec. || introduces 
the Green function in real space and imaginary time. In 
Sec. [I] an analytic solution for the homogeneous electron 
gas is presented. I derive asymptotic formulas for real 
metals and insulators in Sees. IV and The results 



are summarized in Sec. VI together with a discussion. 
Atomic units are used throughout. 



and analogous for the electron part. Here R labels the 
lattice vectors, z«n(r — R) indicates the Wannier orbitals 
corresponding to a set of Bloch wave functions '0„k(r), 
and the integrals extend over the first Brillouin zone with 
volume Q.B- Wannier orbitals are widely accepted to ex- 
hibit exponential localization. In addition to a rigorous 
proof for isolated bands,E3 there is compelling numerical 
eviderijGe to this effect for the general case of composite 
bandsEil Crucially, the localization persists even if. de- 
fects or surfaces break the periodicity of the solid. lj As 
the decay length of Wannier orbitals is of the order of 
the interatomic spacing, only a small number of terms 
contributes significantly to the sum (||) . 



II. THE GREEN FUNCTION 

The kaaginary-time or thermal Green function is de- 
fined aaH3 

g{vy-T) = e(-r)(e(^-''*)"^t(i.')e-(^"'^^)-^(r)) (1) 

where V'^(r) and ijj{Y) denote the creation and annihila- 
tion operator for an electron at r. The Hamiltonian H 
is modified by the chemical potential and the electron 
number operator iV, the angular brackets indicate the 
thermodynamic average, and 6(t) denotes Heaviside's 
step function. Spin degrees of freedom are suppressed to 
simplify the notation. At finite temperature the real vari- 
able r is restricted to — /? < t < with (3 = I/UbT for 
the hole part and < t < /3 for the electron part, and the 
Green function is periodically repeated. Due to the addi- 
tional antiperiodicity Q[t) = —Q{t + I3) it suffices to ex- 
amine the hole part. At zero temperature the period be- 
comes infinite, however, and the properties of both parts 
must be studied separately. Note that Refs. employ 
an alternative prefactor convention G(ir) iQ{—T) for 
the imaginary-time Green function at zero temperature, 
which stems from analytic continuation in the time do- 
main rather than the frequency domain. 

For real materials it is convenient to rewrite the hole 
part of the Green function in the form 



a(r,r';T) =^e„(r,r';r), 



(2) 



a„(r, r'; r) = ^ u;„(r - R)i^„(R - R'; t)wI{v' - R') 

R,R' 



Wn{v - R) = /e-^''-^Vnk(r) 

^l-B J 



(3) 



(4) 



III. THE HOMOGENEOUS ELECTRON GAS 

For the homogeneous electron gas, the hole part of the 
Green function is given by 



Q{r, t) 



1 



g-(efc-Ai)T 
1 -)_ (,(^k-ti)P 



(6) 



with efc = ■ The range of integration can be extended 
to infinity, because the spherical Fermi surface is fully 
contained inside the Brillouin zone. The angular integrals 
are straightforward and lead to 



Q(r, t) 



ke 



ikr 



1 -I- g(efc-M)/3 



dk. 



(7) 



The remaining integral is most simply solved by closing 
the contour across the upper complex half-plane. The 
integrand has relevant first-order poles at ±fc; -f- iji with 



7i 



(8) 



(9) 



and the fermion Matsubara frequencies uji = {21 + 1)tt/ (3. 
Evaluation of the residues yields the final result 



^ oo 

Q{r, t) = 3- e"''"" cos(fc;r - w/t), 

^ 1=0 



(10) 



which explicitly shows the periodicity along the imagi- 
nary time axis. At finite temperature the Green function 
is essentially a superposition of exponentially damped os- 
cillations, whose long-range behavior is dominated by 
the I = contribution. For sufficiently small T the 
chemical potential approximately equals the Fermi en- 
ergy ep = ^^F- The damping constants are given by 
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7/ ~ uji/kp in this regime, and the wave vectors of the 
corresponding osciUations become ki ~ kp. The overall 
decay constant is hence proportional to the temperature 
and given by 7 = 70 = irkBT/kp. 

As the limit T is nontrivial, the Green function at 
zero temperature is most simply calculated from Eq. (|^) 
with a step function instead of the Fermi distribution. 
The integral can be solved analytically and yields 



G{r,T) 



sm^kpr) 



1 



27r2rT 2(27r|r|)3/2 

„ / kpT + ir\ 
ern , 




For the electron part erfc(z) = 1— crf(z) replaces erfi(2:) = 
eTi{iz)/i. The limiting formulas are identical in either 
case and follow from the asymptotic representation of the 
error function erf(z) in the complex plane. For r — » 00 
the Green function becomes 



C/(r,T) 



kp cos{kFr) 



27r2r2 



L-3fsm{kpr), (12) 



which decays algebraically rather than exponentially. On 
the imaginary time axis the Green function decays like 



, , sinikpr) cos(kpr) 
y[r,T) — —- — - 



(13) 



for large positive or negative r. 



IV. METALS 

The explicit solution for the homogeneous electron gas 
illustrates more general principles that also apply to real 
metals. Bands above or below the Fermi level yield con- 
tributions to the Green function analogous to the case 
of insulators, which is discussed below. They are short- 
ranged, so that for |r — r'| 00 the sum (^) is dominated 
by terms arising from partially occupied bands that cross 
the Fermi level. Furthermore, the strong localization of 
the Wannier orbitals implies that the spatial decay prop- 
erties of the relevant Qn are in turn determined by Fn- If 
the homogeneous electron gas is used for guidance, then 
the decay length I/7 for sodium with an average density 
rg = 3.99 (Ref. ^6]) at room temperature T = 298 K may 
be estimated as 85.8 A. In contrast, the decay length of 
the temperature-independent Wannier orbitals is com- 
parable to the interatomic spacing 3.72 A in the bcc cell 
and more than an order of magnitude smaller. The high- 
temperature region where Fn becomes sufficiently short- 
ranged such as not to dominate the asymptotic behavior 
is irrelevant for solid-state physics, because it lies above 
the melting point of the crystal lattice. 

Fn has the mathematical form of a three-dimensional 
Fourier transform and is long-ranged because of rapid 
oscillations of the integrand in Eq. (||). Quantitatively, 



at finite temperature the integrand falls from near unity 
to zero in a narrow region around the Fermi surface. Its 
width is given by 5k cx A:BT/|Ve„k|, where the gradient 
is taken at the Fermi level. Directional effects due to 
the anisotropic crystal structure are ignored here, and 
it is assumed that the energy dispersion is analytic and 
smooth on the scale of 5k. For sufficiently small T the 
latter condition is. always fulfilled. According to basic 
Fourier analysis,EZI the extent of this rapid variation is 
directly proportional to the exponential decay constant 
of the transform, so that 7 oc (Sfc for Fn and hence for Q. 

At r = the integrand changes qualitatively, because 
the Fermi distribution has a discontinuity at the Fermi 
surface. As the Fourier transform of a discontinuous 
function decays algebraically rather than exponentially,LZI 
both the electron and the hole part of the Green function 
now assume the asymptotic form t/(r,r';T) ^ |r — r'|^'' 
with 7y > 0. Their decay on the imaginary time axis is 
likewise algebraic Q{y,t']t) ~ \t\~'^ with k > 0, because 
the exponential exp(— enkT") in the integrand reaches the 
Fermi level and thus exhausts the prefactor exp(/ir). 



V. INSULATORS 

Temperature effects are negligible for insulators if the 
thermal energy ksT is much smaller than the band gap 
A. As this condition is almost always fulfilled below the 
melting point, it suffices to examine the case T = 0. If 
no band crosses the chemical potential, then all Fn are 
spatially short-ranged. For r the integral ^ has 
the straightforward solution F„(R — R';0) = (Jr.r/. At 
finite imaginary times the only significant contribution 
comes from the vicinity of the band edges. The energy 
dispersion can hence be replaced by its harmonic approx- 
imation |e„k - /^l ~ |A -I- i ^^^i,] f^ii'"^*)ij^kj, where m* 
denotes the effective mass tensor at the top of the va- 
lence band for the hole part and at the bottom of the 
conduction band for the electron part. The range of in- 
tegration may now be extended to infinity, and Fn is 
readily seen to be of Gaussian type. The overexponential 
falloff implies that for insulators the asymptotic proper- 
ties of the Green function in real space are determined 
by the Wannier orbitals, which are exponentially local- 
ized. In the weak-binding case, where the band struc- 
ture can be obtained by perturbation from that of the 
homogeneous electron gas, their decay rate is given by 
7 cx aA, where a is the lattice constant.El Many com- 
mon semiconductors have gaps substantially smaller than 
their band width and are thus expected to fall into the 
weak-binding regime, e.g., for silicon A — 1.15 eV (Ref. 
^ ) as opposed to ^(27r/a)2 = 5.1 eV. The proportion- 
ality factor is direction-dependent but about unity, so 
that the decay length I/7 may be estimated as 1.22 A, 
which is of the order of the interatomic spacing 2.35 A 
in silicon, as-expected, and compatible with numerical 
calculations .L3 In the tight-binding regime the Wannier 
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orbitals also decay exponentially, but the dependence of 
7 on the gap is indetefiminate and subject to details of 
the atomic potential.Ell This case may apply for large- 
gap Mott-Hubbard insulators like NiO but is less rele- 
vant in practice, because the GW approximation any- 
way shows sepious deficiencies for such strongly corre- 
lated systems.Ea The falloff on the imaginary time axis is 
exponential with a decay constant 6 = both for the 
electron and the hole part, because the finite difference 
between the chemical potential and the band edges im- 
plies an incomplete cancellation of the prefactor exp(/xr). 

VI. SUMMARY 

In this paper I have investigated the asymptotic prop- 
erties of the one-particle Green function in real space 
and imaginary time. An analytic solution for the ho- 
mogeneous electron gas shows the Green function as a 
superposition of damped harmonic oscillations with a 
wave vector approaching fcj? as T — > 0. The damping 
Q exp(— 7|r — r'|) is exponential at finite temperature 
but changes qualitatively to algebraic falloff ^ |r— r'|^'' 
at T = 0. While the thermal Green function is by defi- 
nition periodic on the imaginary time axis with the pe- 
riod 2/3, it decays like G ~ kl"" when the periodicity 
is lifted at T = 0. I have shown that these results also 
apply to real metals. In contrast, for insulators at zero 
temperature the Green function decays exponentially like 
Q ~ exp(— 7|r — r'|) in real space and Q ~ exp(— 6|t|) in 
imaginary time. 

The limiting formulas can be exploited in electronic- 
structure calculations based on a space-time representa- 
tion and allow an analytic treatment of the long-range 
tails. The asymptotic behavior of derived propagators 
within a perturbation scheme follows directly from the 
Green function. In particular, the self-energy in the GW 
approximation is the product of the Green function Q 
and the screened Coulomb interaction W. In insula- 
tors screening is incomplete with W ~ l/e|r — r'| at 
large distances, where e is the dielectric constant, but 
the self-energy still decays exponentially in real space 
because the Green function does. On the other hand. 
for the homogeneous electron gas at zero temperatureE3 
W ~ cos(2A:i?|r— r'|)/|r— r'l^, so that the self-energy only 
decays algebraically. This distinct asymptotic behavior 
should be appreciated if the electron-gas self-energy is ap- 
plied to real materials. The problem is naturally solved 
by the ad hoc inlxoduction of a gap in the otherwise 
metallic spectrum.E^I 

Finally, the generic dependence 7 oc fcBT/|Ve„k| and 
b — ^ksT for metals as well as 7 oc aA and b = 
for insulators can be used to assess the scaling of dif- 
ferent algorithms for self-energy calculations in the GW 
approximation. If all other parameters remain fixed, then 
conventional reciprocal-space implementations scale like 
N^/T^ for metals at finite temperature and iV^/A^ for 



insulators at zero temperature, whereas the space-time 
method scales like N"^ /T^ and iV^/A^, respectively. This 
difference in efficiency gain as the temperature or the 
band gap increase is an important factor to be taken into 
account for benchmarking purposes. 
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